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o ■ 

^vj . The chargino/neutralino pair production is one of the benchmarking processes of ILC. 

These processes are interesting not only because it allows high precision measurement 

^ , of chargino and neutralino masses, but also for the reason that the separation of W 

H— 5 . and Z bosons through their hadronic decay products requires excellent jet resolution 

ly-N ' being a good benchmark of the detector performance. The analysis based on the SiD 

detector concept with four jets and missing energy final state will be presented. The 

uncertainty of chargino and neutralino cross sections can be determined with precision 

QJ | of 0.9% and 4.2% respectively. The mass uncertainties are obtained with a template 

r \^ . fitting method achieving precision of better than 1 GeV. 

1 Introduction 

q ■ The chargino/neutralino pair production or so-called 

• T* \ SUSY Point 5 is one of the benchmarking processes 
^►y proposed for the ILC pQ. For the SUSY parame- 

,J^ \ ters listed in Table [TJ the lightest chargino mass is 

i Mh - smaller than slepton masses but at the same time is 

much larger than the lightest neutralino mass. Xi is 

i-H ■ the Lightest Supersymmetry Particle (LSP), there- 

fore charginos dominantly decay into LSP and on- 

~*~ ■ shell W bosons. In this scenario the x® mass is almost 

^L \ degenerate with the chargino mass and it similarly 

{-*\ ■ decays into LSP and Z bosons. 

1^. ■ _i_ _ ~j ^o-Ott^+tt^- Table 1: SUSY point 5 parameters. 

^ ■ e+e -> xi Xi -► XiXW^W 

e + e- -► xlxl "+ X?X?Z°Z° 
^H . 

If we consider all-hadronic decays of the gauge bosons in the final state, the above 

• tH . processes will both be defined by a signature of four jets with large missing energy. Once 
/\ ' the W or Z bosons are successfully reconstructed, their energy distribution will provide 

information on the chargino and neutralino masses. The sensitivity of the mass measurement 
is determined by template fitting. Since this measurement requires separation of pairs of 
W and Z bosons reconstructed from two jets, fine jet resolution and a good performance of 
Particle Flow Algorithm (PFA) [2 are essential. 

A Monte-Carlo sample with y/s of 500 GeV and integrated luminosity of 500 fb _1 was 
used for this study. It consisted of the chargino/neutralino signals and SUSY backgrounds 
such as e + e~ — >• X1X2 an d slepton pair production. Apart from this sample, a few more 
samples with the same statistics were generated in order to estimate the uncertainty of the 
mass measurement, each with one SUSY particle mass shifted by certain amount (Table EJ). 
A Standard Model background sample at y/s = 500GeV was generated separately. All the 
samples have the same polarization, namely right handed 80% electrons and left handed 30% 
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parameter 


value 


m 


206 GeV 


m 1/2 


293 GeV 


tan /? 


10 


A 





/i 


375 GeV 


M xl 


115.7 GeV 


M xt 


216.5 GeV 


M x% 


216.7 GeV 



positions. The cross-section of chargino production with all hadronic final state is 57.1 fb 
and the neutralino one is 11.0/6. Note that the cross-section of neutralino pair production 
is much smaller than the chargino one, therefore the separation of the two, especially the 
neutralino selection in the presence of a large chargino background is challenging. 



Template M^o (GeV) M~± (GeV) M^o (GeV) 



standard 


115.7 


216.5 


216.7 


X? + AM 


115.7 + AM 


216.5 


216.7 


Xi +AM 


115.7 


216.5 + AM 


216.7 


X^ + AM 


115.7 


216.5 


216.7 + AM 



Table 2: MC sample parameters. AM = -0.5/0.5/2 GeV or 6 GeV (only for xl) 



2 Analysis Methods 

All samples used in this study were generated using the Monte Carlo program Whizard 
1.4 |3]. PYTHIA 6.4 [4] is used for parton showering, fragmentation and decay to provide 
final state observable particles. The detector response to the generated events is simulated 
by the Linear Collider specific code SLIC [5] based on the Geant4 toolkit [7] [8]. The full 
detector description of the baseline Silicon Detector (sid02 [9!) is employed. 

The particle reconstruction was performed within the lcsim 1.5 framework [6]. Pattern 
recognition on hadronic final states used a PFA algorithm devised for SiD pj]. The jet clus- 
tering and further physics analysis were performed within the MarlinReco framework [TO] . 
The Durham algorithm [12] was used for jet finding and the jet number was forced to be 
four rather than setting a y cut value. 

The analysis strategy is defined as follows. Common preselections for the xf x^~ an d X2X2 
events were defined to suppress the Standard Model background, as described in Section I2TT1 
The cross-sections of both signal events could be determined after they were separated using 
the dijet mass correlation as described in Section [2^21 The energy distributions of the gauge 
bosons as decay products gave information on the SUSY particle masses. Uncertainties of 
the masses from this measurement were determined by the template fitting method presented 
in Section SI 

2.1 Signal Preselection 

The analysis started with selection of all-hadronic channel events. All reconstructed particles 
were clustered into four jets. A series of selections were applied (Table [3]). The selection 
criteria are explained briefly as following: 

• Total visible energy: Events with missing energy are selected; 

• Number of tracks: Leptons reconstructed as jets with low multiplicity are eliminated; 

• Thrust; cos9 t h rust : Jets originated from signal events should be distributed more 
uniformly; 

• Ej e t] fraction of EM energy in each jet; lepton energy in each jet: Ruling out jets 
falsely reconstructed from energetic leptons; 
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• 0(i,j): opening angle between jet i and j (ordered by energy); 

• Acoplanarity: The WW or ZZ pair from XiXi or X2X2 decay are acoplanar, different 
from WW pairs in SM. 

cut value 

Jet number = 4 

Total visible energy < 250 GeV 

Number of tracks > 20 

Thrust < 0.85 

COS Othrust < 0.9 

E jet < 10 GeV 

Fraction of EM energy in each jet < 80% 

(9(1, 2) > 60° 

0(1, 3), 0(1, 4) , 61(1, 3) > 40° 

(9(2, 4), 0(3, 4) > 20° 

lepton energy in jet 1 < 40 GeV 

lepton energy in jet 2 < 40 GeV 

lepton energy in jet 3 < 30 GeV 

lepton energy in jet 4 < 20 GeV 
Acoplanarity of two reconstructed gauge bosons > 10° 

Table 3: Pre-selections. 0(i,j) defines the opening angle between the ith. and jth jet. Jets 
were ordered in energy, e.g. jet 1 is the most energetic jet. 

As a result the Standard Model background was well suppressed (Figured]). 

2.2 Chargino/Neutralino Signal Separation 

To reconstruct two gauge bosons using the four jets we have to determine the jet pairing 
corresponding to correct combinations. It is done by minimizing the quantity of (m(ji, J2) — 
m w/z) 2 + ( m (J3? Ja) — mw/z) 2 , where m(ji, j'2) is the invariant mass of the jets j\ and J2- 
All jets were ordered by energy and all possible jet permutations were considered. The W 
mass was used when selecting XiXi events, and the Z mass for X2X2 events. The correlation 
of the two reconstructed boson masses is shown in Figure O Reconstructed boson masses 
populate different regions in the histogram for chargino and neutralino events, and their 
separation can be achieved by a selection indicated by a thick line in Figure [2j Events in 
the lower left part (mwi + mw2 < 172 GeV) were classified as XiXi events and in the 
upper right part (mwi + mw2 > 172 GeV) as X2X2 events. For chargino events selection 
we also add a selection to remove events in the most lower left corner of 2D mass histogram 
(mwi + wtW2 > 130 GeV^), which are most likely due to the remaining leptonic or semi- 
leptonic decays of the bosons. 

2.3 Kinematic Fitting 

In the rest frame of X\ /X2? W^/Z bosons are monochromatic and their energy is determined 
by three masses: m~± /m^o, m w ±/mz and m^o. For a boosted parent Xi /X25 ^ ne boson 
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Figure 1: Examples of SM background suppression cuts for the chargino selection. Solid 
lines are the chargino signal and dashed lines are the Standard Model background. Left: 
Total visible energy; Right: Acoplanarity of two reconstructed W bosons. 
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Figure 2: Reconstructed boson masses from the four jets, for chargino selection. Left: pure 
chargino signal; Right: pure neutralino signal. The region between the two straight lines 
shows an allowed chargino selection window. 
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Figure 3: Mass (left) and energy (right) distributions of reconstructed W boson before 
(black) and after (red) kinematic fitting. 



energy distribution depends on the SUSY particle masses which will determine the edges 
of the energy spectrum. Our approach to obtain these masses is to compare the energy 
spectrum with a reference MC sample template (Table [2j). Therefore the best possible 
energy resolution is important for a high-precision mass measurement. Using the fact that 
the two bosons reconstructed from four jets have the same mass in both XiXi an d X2X2 
cases, kinematic fitting with one constraint (rribosoni = mboson2) can be used for improving 
the energy distribution assuming that the detector resolution is known. Kinfit in Marlinreco 
package [13] is used for the fitting. The comparison of reconstructed W mass with and 
without fitting is shown in Figure 03 

3 Cross-section Results 

The reconstructed boson energy distribution are shown in Figure 2] after preselections and 
signal separtion. The chargino sample contains 15362 signal events, 1997 events of neu- 
tralino pair background, 65 events of other SUSY background and 2980 events of Standard 
Model background. Similarly the neutralino sample contains 1659 events of pure neutralino 
signal, 2395 events of chargino background, 9 events of other SUSY background and 865 
events of SM background. The uncertainty of the cross section measurement for chargino 
and neutralino pair production with all hadronic final state is 0.9% and 4.2% respectively 
neglecting the uncertainty of efficiencies. 

4 Mass Uncertainty and Template Fitting 

The determination of mass sensitivity can be explained on the example of the Xi^XT process. 
The W energy distribution is shown in Figure [5] for different Xi masses. The were all 
normalized to the same cross section to explore the pure kinematic effect. The change of 
upper and lower edges of the spectrum after normalization is clearly visible. Assuming that 
we can generate a lot of templates, the mass of the SUSY particle masses can be obtained 
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Figure 4: The inclusive histogram of reconstructed W energy for XiXi selection (left) and 
Z energy for X2X2 selection (right). 



by comparing the boson energy spectrum from the data with those of various templates. 
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Figure 5: Left: W energy distribution for different Xi masses; Middle: W energy distribution 
normalized to the same cross-section; Right: Difference of the normalized W energy spectrum 
with respect to the standard template (M^o = 115. 7 GeV). 

In this analysis, one 'data' sample is taken from the stan- 
dard template (Tabled]) without any SM background. The x 2 
of a template with respect to the 'data' is defined as following: 



Nbins / ^2 

2 _ V~^ {JJ tempi at e,i Vdata,i \ u ) 



£ 



(i) 



-q template,! _r data,i 



J SM,i 



where t/...^ denotes the content of the ith bin in the boson 
energy histogram. The S term is added as a Gaussian smear- 
ing of the central value of y te mpiate,i ~ Vdata,i to take the SM 
background into account. The x 2 dependence on the mass 
is parabolic, with the minimum located at the actual mass. 
The half- width of the x 2 = X m i n + 1 gives one a uncertainty 
on the mass [14] as illustrated in Figure [6] The statistics of 



X min + 




Figure 6: Sketch of tem- 
plate fitting method (not to 
scale). 
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the templates are much larger (~ 10 times) than the 'data'. 
The uncertainty results are stable for different binning of the 
histogram, and the minimum Xmin over degree of freedom should be close to one. 

Figure [3 shows the x 2 parabola from the W energy spectrum, varying xf an d Xi masses 
respectively. It should be noticed that the point for AM = is calculated for an independent 
standard template other than the 'data' sample. Similarly for the X2X2 events the x 2 
parabola can be obtained from various X2 an d X? masses (Figure [8)). AM has the the value 
of -0.5, 0, 0.5 and 2 GeV for all target masses, and 0.6 GeV only for x°. This is because the 
X2X2 sample has a worse purity comparing to xf Xi due to a large background of the latter, 
and therefore the change in Z energy edge is less sensitive to the X2 m &ss. These figures are 
obtained with the energy bin width of 5 GeV. 
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Figure 7: The x dependence on the shifted mass of xf (left) and Xi (right), calculated for 
the xt Xi signal. 

I Ne2 Mass Fit (Neu2Neu2 events) rebin 5~| I Nel Mass Fit (Neu2Neu2 events) | 
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Figure 8: The x dependence on the shifted mass of X2 (left) and Xi (right), calculated for 
the X2X2 signal. 

Uncertainties of SUSY particle masses can be estimated as in Table H] from the parabolas 
above. These are taken after averaging from different binning, namely energy histograms 
with 1 GeV, 2 GeV, 5 GeV or 10 GeV per bin. 

Each uncertainty was calculated by varying that specific mass alone, assuming other 
SUSY particle masses are measured perfectly. In order to estimate possible correlations 
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X^Xi 


xlxl 


AM(xf ) 472 MeV 
AM(xi) 156 MeV 


AM(x°) 975 GeV 
AM(x?) 279 MeV 



Table 4: Mass uncertainty results. 

between the mass uncertainties, several templates are needed with masses, e.g. xt/Xi or 
X2/X1 shifted simultaneously. More samples will be generated for further study. 

5 Conclusion 

The chargino/neutralino pair production at the ILC with full SiD detector simulation are 
studied. The cross-section of XiXi an d X2X2 are measured with sensitivity to 0.9% and 
4.2%. Template fitting method is applied to determine the uncertainty of the chargino and 
neutralino masses, and the mass uncertainties obtained were all below 1 GeV. 
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